Classical Trajectory Perspective on Double Ionization Dynamics of Diatomic 
Molecules Irradiated by Ultrashort Intense Laser Pulses 
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In the present paper, we develop a semiclassical quasi-static model accounting for molecular dou- 
ble ionization in an intense laser pulse. With this model, we achieve insight into the dynamics of two 
highly-correlated valence electrons under the combined influence of a two-center Coulomb potential 
and an intense laser field, and reveal the significant influence of molecular alignment on the ratio 
of double over single ion yield. Analysis on the classical trajectories unveils sub-cycle dynamics of 
the molecular double ionization. Many interesting features, such as the accumulation of emitted 
electrons in the first and third quadrants of parallel momentum plane, the regular pattern of corre- 
lated momentum with respect to the time delay between closest collision and ionization moment, are 
revealed and successfully explained by back analyzing the classical trajectories. Quantitative agree- 
ment with experimental data over a wide range of laser intensities from tunneling to over-the-barrier 
regime is presented. 

PACS numbers: 33.80.Rv, 34.80.Gs, 42.50.Hz 
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I. INTRODUCTION 

Atoms and molecules exposed to high-intensity ultra- 
short laser pulses have been attracting much attention 
during the past ten years and plenty of striking phe- 
nomena emerged such as multiphoton ionization (MPI), 
above threshold ionization (ATI), high order harmonic 
generation (HHG) and nonsequential double ionization 
(NSDI) [J, [|. The dynamics behind the NSDI is most 
complicated because it deals with two highly entangled 
valence electrons Q. In contrary to the situation of sin- 
gle ionization where the experimental data can be ac- 
curately predicted by the single active electron (SAE) 
model [4ij , the NSDI data are usually higher by many or- 
ders of magnitude than the sequential Ammosov-Delone- 
Krainov (ADK) tunneling theory that bases on the SAE 
approximation. This fact suggests that the electron- 
electron correlation plays a crucial role in NSDI physics. 
Diverse correlation mechanisms like the shake-off [5| and 
the recollision were proposed to account for the elec- 
tron correlation in the experimental observations. It is 
commonly believed now that the rescattering is the dom- 
inant mechanism for the excessive atomic double ioniza- 
tion (DI) yields 0. 

Compared to atomic case, the NSDI in the molecules is 
more complicated because of the additional freedom that 
molecules have0, E3|- In addition to the typical strong- 
field phenomena like excessive DI yields in the tunnel- 
ing regime, high harmonic generation of the driving laser 
field, and momentum correlation between the two emit- 
ted electrons, some new phenomena such as bond soft- 
ening [ll[, zero photon dissociation [l2j], and alignment 
dependence of DI yield [H, EH EH were observed in re- 
cent molecular experiments. 



However, these experimental data for molecular dou- 
ble ionization is far from well understood in theory. The 
complex dynamics of correlated e — e pair responding to 
two-center nuclear attraction and laser force poses a great 
challenge to any quantum theoretical treatment. For in- 
stance, a fully-dimensional quantum-mechanical compu- 
tation from first principle is very very time expensive 
even for the simpler case of highly symmetric atoms [l6j . 
This leaves approximate approaches developed recently, 
such as one-dimensional quantum model [l7| , many-body 
S-matrix [l8| and simplified classical methods [l9[ . In ei- 
ther case, the complex electron dynamics which is crucial 
for molecular DI is not fully explored and the theoretical 
results can not account for experimental data quantita- 
tively. 

In this paper, we explore a feasible semiclassical the- 
ory, following the treatment in our recent Letter (20j , ex- 
tending the calculations and presenting a more detailed 
account of the theoretical methodology. 

Our calculation is capable of quantitatively reproduc- 
ing the unusual excess DI yield of nitrogen molecules for 
a wide range of laser intensities from 5 x 10 13 W/cm 2 to 
1 x 10 15 W/cm 2 . The significant influence of molecular 
alignment on DI yield is virtually revealed: i) The ratio 
between double and single ionization yield is less for per- 
pendicular molecules than that of parallel molecules; ii) 
This anisotropy becomes more dramatic if the molecules 
are irradiated by a shorter laser pulse. In particular, 
our model provides an intuitive way of understanding 
the complex dynamics involved in the molecular DI with 
back analyzing the classical trajectories. Their sub-cycle 
energy evolution verify that collisions between electrons 
determines the fate of molecular DI, therefore consolidate 
the classical rescattering view of molecular DI. Statics 
based on plenty of individual trajectories show that elec- 
trons are most likely to be emitted at 30° off laser peak, 
indicating the accumulation of emitted electron pairs at 
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±0.5a.u. in the first and third quadrants 



of parallel momentum plane (k± , k}} ) , which is consistent 
with the experimental results. Moreover, the ejected elec- 
trons emerge in the same momentum hemisphere and op- 
posite momentum hemisphere alternately, depending on 
weather the delayed time between closest collision and 
ionization is odd or even half laser cycles. 

Our paper is organized as follows. In Sec. II we present 
our semiclassical quasi-static model. Calculations on 
double ionization yield for nitrogen molecules are pre- 
sented and compared with experimental data in Sec. III. 
In Sec. IV, we give an analysis on the classical trajectories 
and unveil sub-cycle dynamics of the molecular double 
ionization Sec.V is our conclusions and discussions. 



II. MODEL 

We consider a diatomic molecule with two valence elec- 
trons irradiated by a laser field whose temporal and spa- 
tial distribution (see Figflja)) is expressed as 

irt 

e{t) = e (Rl,Zl) sin 2 ( — ) cos(ujt)e z . (1) 
nl 

In experiments, molecules are driven through the laser 
beam perpendicularly. Thus the external field along the 
propagation direction of the laser beam is approximately 
constant and in the lateral direction it can be treated 
as an ideal Gaussian beam, i.e., eq(Rl, Zl) = eq(Rl) = 
eq exp (—R l /Rq), where £o is the peak laser field, Rl 
represents the position of the molecules in the laser beam, 
and Rq is the radius of the beam. The maximum of Rl 
is chosen to be three times or more of Rq to ensure the 



convergence of our results. A 



enveloped laser with 



the full width of nT is used in our calculations, where T 
and w are period and angular frequency of the laser field 
respectively and n denotes the number of optical cycles. 

In quantum mechanics, ionization processes are de- 
scribed by the evolution of wave packets, however, in our 
model they are traced by launching a set of trajectories 
with different initial parameters instead (To some extent, 
the distribution of these parameters can be regarded as 
an alternative emulation of quantum mechanics). For a 
certain trajectory, the evolution of two electrons can be 
described by classical Newton equations: 



dh 



dt 2 



£(i)-VOC + K e ). 



(2) 



Here e(t) is the external laser field discussed above. The 
index i denotes the two different electrons. V* e and V ee 
are Coulomb interaction between nuclei and electrons 
and between two electrons, respectively. 
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FIG. 1: (color online), (a) The spacial and temporal profile 
of a typical laser beam, (b) The coordinate system adopted 
in our paper. 



where r a i and rb% are distances between the ith electron 
and nucleus a and b, respectively (as shown in Fig[T](b)). 

To solve the Newton equations, we need the initial con- 
ditions of the two electrons, including the initial time, ini- 
tial position and momentum. Their distribution are quite 
different for the tunneling and over-the-barrier regime 
due to the different ionized mechanism of the first elec- 
tron. 



A. Tunneling regime 

In the long-wavelength limit, the laser field varies 
slowly in time and can be regarded as a quasi-static field 
compared with valence electron's circular motion around 
nuclei. Under this field, the Coulomb potential between 
nuclei and electrons is significantly distorted. When the 
instantaneous field (at time to) is sufficient strong but 
still smaller than a threshold value (see FigJ^a)), one 
electron is released at the outer edge of the suppressed 
Coulomb potential through quantum tunneling with a 
rate w(to) given by molecular ADK formula [21j. 

The electron tunnels out through a saddle point (l9l| 
directing to a channel of the local minimum in the com- 
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FIG. 2: (color online), (a) Tunneling ionization, (b) Over- 
the-barrier ionization, (c) The contour plot of the combined 
Coulomb potential and external laser field. It clearly shows 
that the saddle point locates approximately along the direc- 
tion of the external field. 



bined potential of the Coulomb interaction and the exter- 
nal field (see, Fig[2Ic)). Because the difference between 
the direction of the saddle point and the external field 
is very small, we safely regard the external field direc- 
tion (z axis) as the tunnelling direction. Thus, the initial 
position of the tunnelled electron can be derived from 
following equation, 



1 



1 



*(r) 



-n-dr + I pl ~ zis(t ) = 0, (4) 



with x\ = yi = 0. The wavefunction 'F is given by the lin- 
ear combination of the atomic orbital-molecular orbital 
(LCAO-MO) approximation. Taking iV^~ for example, we 



choose 



A 3/2 



-\r 



as the trial function to construct 



the molecular orbital 'F(r) = c[^>(r 2)+(^(r(,2)], where c is 
the normalization factor. The parameter A, which equals 
to 1.54 for N2 , is determined through variational ap- 
proach. That is, calculate the variational energy for the 
given wavefunction and assume it equals to the second 
ionization energy of the molecule. The initial velocity of 
tunnelled electron is set to be (v± cos ip, v± sin tp, 0), with 
v± having the same distribution as that in atomic case, 
i.e., 



w(v±)dv± 



2{2I P i) 1/2 



e(to) 



v± v 
exp(- 



2 m P i) 1/2 



e(*o) 



■)dv ± , (5) 



where tp is the polar angle of the transverse velocity uni- 
formly distributed in the interval [0, 2ir}. 

In order to get the initial velocity distribution, we 
employ a technique widely used in classical Monte 
Carlo(MC) simulation. Firstly, we generate two ran- 
and Wtest in the interval [0, u™ ax ] and 



dom number u^ st 



[0, Wmax], respectively. If w(^ esi ) > w test , 



is kept 



as the initial transverse velocity, otherwise it is rejected 
and the above procedure is repeated. 

For the bound electron, the initial position and mo- 
mentum are depicted by single-electron microcanonical 
distribution (SMD) 



F(r 2 , P2 ) = kS[I p2 - p 2 2 /2 - W(r a2 ,r b2 )}, 



(6) 



where k is the normalization factor, I p2 denotes the 
ionization energy of molecular ions such as , and 
W{r a2 ,r b2 ) = —l/r a 2 — 1/V&2 is the total interaction po- 
tential between the bound electron and two nuclei. 



B. Over-the-barrier regime 

The above scheme is in the spirit of our quasi-static 
model for atomic DI [24| and is applicable only when 
the instantaneous laser field is weaker than the threshold 
value [25j . To give a complete description of the DI of 
molecular system for the whole range of laser intensities 
(see FigH]), one needs to further extend the above model 
to over-the-barrier regime (FiglJb, in such case Eq.(j4| has 
no real roots). This is done by constructing the initial 
conditions with double-electron microcanonical distribu- 
tion (DMD) [H, i.e., 



F(ri, ,r 2 ,pi,p 2 ) 



1 



[/a(ri,Pi)//3(r 2 ,P 2 ) 



//3(ri,Pi)/a(r 2 ,P 2 )], 



(7) 



with 



/a,/3(r,p) = k S[I pl 



2 



W(r a ,r b )-V a ,p(r)}, (8) 



where V^.^r) represents the mean interaction between 
two electrons, V a ^(r) [1 — (1 + Kr b ^ a )e~ 2Krb - a ]. k can 
be obtained by a variational calculation of the ionization 
energy of molecules (re — 1.14 for N 2 ). Details can be 
found in Ref.^. 

We would like to give some remarks here, (i) In our 
calculations, the trajectories obtained from DMD are 
weighed by w(to) [21(. (ii) Part of the electrons ob- 
tained from DMD could "self-ionize" even without the 
presence of external field. To avoid such a unphysical 
sclf-ionization, we evolute the electrons for several opti- 
cal cycles freely and abandon those samples whose energy 
are greater than zero during the free evolution process, 
(iii) The initial total energy distribution of the two elec- 
trons is plotted in Fig(3J where we see a long tail on both 
sides. They are cut off by introducing two parameters 
E m i n and E max , which satisfy 



i^::; P (E)dE 



E 



M ■ 



(9) 



where Em ~ I p i + I p2 is the most probable energy and 
p(E) is the state density around E. 
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FIG. 3: Total energy distribution of the two electrons. E m i n 
and E max are used to cut off the long tail. Only the electrons 
with a total energy between E m i n and E max is adopted in our 
calculation. 



With the above initial conditions, the Newtonian equa- 
tions are solved using the 4-5th step-adaptive Runge- 
Kutta algorithm and DI events are identified by energy 
criterion. In our calculations, more than 10 5 weighted 
(i.e. by rate w(to)) classical trajectories of electron pair 
arc traced and a few thousands or more of DI events 
are collected for statistics. Convergency of the numeri- 
cal results is further tested by increasing the number of 
launched trajectories twice. 



III. CALCULATION ON THE DI OF N 2 : 
COMPARISON WITH EXPERIMENTS 
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FIG. 4: Comparison between DI dataflb]] and theory for ni- 
trogen molecule. 0.185 PW/cm 2 is the threshold intensity 
separates the tunneling regime and over-the-barrier regime 
as schematically plotted in FigO To our knowledge, the re- 
sults from our theoretical model are the first to be in good 
agreement with experimental data for a wide range of laser 
intensities from tunneling regime to over-the-barrier regime. 
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In this section, we apply our theory to the nitrogen 
molecules and compare our model calculation with ex- 
perimental data. The parameters are chosen to match 
the experiments, e.g. the internuclear separation is 
2. 079a. u., the first and second ionization energy are 
Ipi = 0. 5728a. u., I P 2 — 0. 9989a. u., respectively. The 
laser frequency is u = 0. 05695a. u. corresponding to a 
wavelength of 800nm. Based on our model, we success- 
fully reproduce the plateau representing the excessive DI 
yields due to the correlation between electrons , and re- 
veal a slight shallow dip at intensity of 0.4 x l0 15 W/cm 2 
indicating a transition point below which the electron 
correlation plays role. In addition, our model calcula- 
tions confirm the dramatic dependence of DI yield on 
molecular alignment. In FigJH the ratio between dou- 
ble and single ionization yield is plotted with respect 
to the peak laser intensities from 5 x 10 13 W/cm 2 to 
1 x 10 15 W/cm 2 . The number of optical cycle is cho- 
sen as 37 to match the experiments of Cornaggia flpj ] . 
Our numerical results show good agreement with the ex- 
perimental data over the whole range, thus confirm the 
validity of our model. The curve can be clearly divided 



FIG. 5: The molecular alignment dependence of DI ratios for 
laser intensity of 0.15PW/cm 2 . 



into three parts: the low intensity regime (5x 10 13 W/cm 2 
to 1.85 x 10 14 W/cm 2 ), where tunneling ionization plays 
an important role; the very high intensity regime (5 x 
10 14 W/cm 2 to 1 x 10 15 W/cm 2 ) which can be also well 
described by ADK theory, indicating that the electrons 
are pulled out sequentially. In particular, there exists a 
plateau (moderate intensity) regime with a shallow dip, 
which can be seen from the corresponding experiments 
in Ref.[ic| Fig. 10. Back analysis about classical trajecto- 
ries indicates that the underlying dominant mechanism 
in these three regimes are widely divergent as will be 
shown latter. 

In our producing data in FigU]the molecular alignment 
is set to be random and our result is obtained by averag- 
ing over different orientations. However, the inherent nu- 
clear degrees of freedom of molecules do manifest them- 
selves as the significant alignment effect in our model. 
On the other hand, recent progress in experimental tech- 
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nique has make it possible to control molecular alignment 
by applying a weak pre-pulse (la |. With these consider- 
ations and further application of our model, we calculate 
the ratios between double and single ionization accord- 
ing to different angles between molecular alignment and 
laser polarization. Main results are presented in Fig|5] 
It clearly shows two tendencies, i) The ratio between 
double and single ionization yield is less for perpendic- 
ular molecules than that of parallel molecules; ii) This 
anisotropy becomes more dramatic when the molecules 
are irradiated by a shorter laser pulse. 

Compared with shorter pulse case, Figj5]Ja) shows the 
ratios between double and single ionization keeps almost 
a constant until 35° and then decreases rapidly. At 90° 
(i.e., perpendicular case), the ratio become half of that 
in parallel case. This result is qualitatively in agree- 
ment with experimental data (l5j . Quantitatively, our 
factor 2 is larger than the observed ratio of 1.1 between 
parallel and perpendicular case [lol ]. This discrepancy 
comes from the unperfect control of molecular alignment 
in practical experiments due to the rotation of molecules 
and many other uncertain factors. Based on our the- 
oretical results, a simple estimation will conclude that 
only about 15% of the molecules are "switched" to the 
appointed direction additionally. Further explorations 
show that molecular alignment also significantly affects 
the correlated momentum distribution of emitted elec- 
trons. Details will be presented elsewhere [22j]. 



IV. SUB-CYCLE DYNAMICS OF MOLECULAR 

DI 



tunnelling regime, l<0.1 85PW/cm 
1.oRb) 




-0.5 



20 21 



22 23 



CD 
£Z 

LU 



ii. plateau regime, 0.185Pw/cm 2 <l<0.5PW/cm 2 




iii. sequential ionization regime, l>0.5PW/cm 
12 r 




21 22 23 24 

Optical cycle (t/T) 



FIG. 6: (color online). Typical energy evolution of the elec- 
tron pair in three different regimes: (i) tunneling regime, (ii) 
plateau regime and (iii) sequential ionization regime. Verti- 
cal dashed lines indicate the moment when significant collision 
between two electrons emerge. 



With our model we are capable to extract vital in- 
formation by tracing back the history of individual DI 
trajectories[27|. Such classical trajectory perspective 
provides us an intuitive way towards understanding the 
complex dynamics involved in molecular DI. 



A. Typical trajectories responsible for emitting 
electron pairs 

In Fig|3l the vertical dot line indicates the threshold 
value of 0.185 PW/cm 2 , which separates the double ion- 
ization data into two parts, the tunneling regime and 
over-the-barrier regime. When peak laser intensity is be- 
low this value, there exist two dominant processes re- 
sponsible for emitting both electrons, namely, collision- 
ionization(CI) and collision-excitation- ionization(CEI), 
as shown in Figj6]^a) (b) , respectively. The unique feature 
of these two mechanism is that the first electron tunnels 
out through the potential barrier and escapes instead of 
being pulled out directly by the laser field. This tunnel- 
ing mechanism can't be included by any classical theory, 
thus we have not choice but place an electron outside the 
barrier directly. As a result, we usually notice that the 
tunneled electron has a higher initial energy as shown in 



FigJSJa) and (b); while in other trajectories, the two elec- 
trons almost have the same initial energy, a reflection of 
the identity. For CI, the tunnelled electron is driven back 
by the oscillating laser field to collide with the bounded 
electron near its parent ions causing an instant (~ at- 
tosecond) ionization. The condition for occurrence of CI 
is usually hard to reach. It requires the kinetic energy of 
the returned electron exceeds the ionization energy of the 
bound electron and they should move close to each other 
to share the excess kinetic energy. However, as pointed 
out by other numerical studies, the returning trajectories 
can be assumed to be a freely spreading Gaussian wave 
packet whose return width is about 30a.u.[28[. So most 
of the cases, the returned electron just slides apart from 
the bound electron and only contribute a little propor- 
tion of its kinetic energy to the inner one, in this case 
CEI occurs instead. For CEI, the inner electron is firstly 
excited through collision with the returned electron, then 
undergoes a time-delayed (~ a few optical periods) field 
ionization of the excited state. During the time inter- 
val, the returned and escaped again electron oscillates in 
the laser field like a near-free electron while the excited 
but still bounded electron move around the cores until 
accumulate enough energy to cause DI. 
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When the instantaneous laser field is above the thresh- 
old value, over-barrier ionization emerges. In this regime 
we observe more complicated trajectories for DI pro- 
cesses. Except for CI (FigHOc)) and CEI (Fig[f|d)) tra- 
jectories similar to tunneling cases, there are multiple- 
collision trajectories as shown in Fig[6je),(f) as well as 
collisionless trajectory of FigJSJg). In FigJSJe) and (f), 
two valence electrons entangle each other initially , ex- 
perience a multiple-collision and then emitted. These 
two kinds of trajectories exits around the threshold in- 
tensity where the over-barrier ionization start to emerge 
but the external field is not strong enough to pull out the 
electron immediately without the assistant of Coulomb 
repulsion between two electrons. The four types of tra- 
jectories indicated by Fig[6f c-f ) represent the dominant 
processes of DI in the plateau regime from 0.185PW/cm 2 
to 0.5PW/cm 2 , which are much more complicated than 
that of tunneling regime, but still accompanied by once 
or multiple times of collisions between two electrons [2!| . 
However, above 0.5PW/cm 2 , the double ionization is 
dominated by a collisionless sequential ionization whose 
typical trajectory is represented by FigJS]Jg). In this 
regime, our numerical results tend to ADK theory. 



B. Time delay effect on correlated momentum 

A key feature distinguishing classical (or semi- 
classical) model from quantum-mechanical treatment is 
the most conspicuous and least ambiguous implication 
of three time scale during the history of individual DI 
trajectories [30| . i.e., the time of single ionization, the 
time of closest collision and the time of DI. The time 
of single ionization is adopted as the start time of our 
simulation, and determines the weight of each trajectory. 
After single ionization, the electron travel much of the 
time in the intense laser field like a classical object until 
recollide with the bound electron, thus lead to DI even- 
tually. Usually, there is a time lag between the time of 
DI and closest collision. Such time delay determines the 
fate of two electrons, and is also firmly related to the final 
parallel momentum. 

Like the atomic simplified model[24[ neglect- 

ing most of the unimportant factors, such as the spacial 
distribution of the laser beam and temporal envelop of 
the laser pulse, would be more convenient while still give 
clear interpretation to most of the phenomena. So in this 
section, we assume the molecules are exposed to a cosine 
external field and the first electrons tunnelled out in the 
first half optical cycle. The electric field is switched off 
using a cos 2 envelope during two optical cycles at last. As 
we would discuss later, many important information can 
be extracted from the distribution of laser phase at the 
moment of closest collision and ionization [3j [H[ . We 
choose three typical laser intensities, i.e. 0.12PW/cm 2 , 
0.4PW/cm 2 and lPW/cm 2 , representing the tunneling 
regime, plateau regime and sequential ionization regime, 
to analyze in detail. 
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FIG. 7: (color online). DI yield vs laser phase when (a) 
the two electrons become closest; (b) both the electrons are 
ionized, at different laser intensity 0.12PW/cm 2 (triangle), 
0.4PW/cm 2 (circle) and IP W/cm 2 (square), respectively. The 
solid lines represent the laser field. 



FiglDJa) shows the diagram of DI yield versus laser 
phase at the moment of closest collision. In the tunnel- 
ing regime (i.e., 0.12PW/cm 2 ), we note that the collision 
can occur throughout most of the laser cycle and the peak 
emerges slightly before the zeroes of the laser field, con- 
sistent with the prediction of simple-man model [fj and 
recent results from purely classical calculation [27j ■ How- 
ever, for other two cases, the collision between two cor- 
related electron turns to occur mainly at peak laser field. 
This is because ionization mechanism changes at the 
transition to over-the-barrier regime, where both elec- 
trons rotate around nuclei initially and they usually be- 
come closest before one of them is driven away by the 
external field (laser peak) rather than during the recolli- 
sion process (zeroes of the laser field). 

FigEUb) confirms that most DI occurs around the max- 
imum of laser field for both tunneling regime and over- 
the-barrier regime. More interestingly, compared to two 
other cases we observe a peak shift of ~ 30° off the field 
maximum for the tunneling case. It is due to the larger 
fraction of CI trajectories in this regime. With assuming 
that the colliding electrons leave the atom with no signif- 
icant energy and electron-electron momentum exchange 
in final state is negligible (these assumptions have been 
checked by directly tracing trajectories), the parallel mo- 
mentum ki 2 of each electron results exclusively from the 
acceleration in the optical field: k\ 2 = ±2y/U p smujti on 
[32l ]. The above shifted peak indicates the accumulation 
of the emitted electrons at = k% = ±0.5a.u. in the 
first and third quadrants of parallel momentum plane 
(fc| , ko). It is consistent with the experimental data of 
Ref.[ljf (see their Fig.2). 

Usually, there is a time delay between closest collision 
and ionization. Fig[8] shows the profound effects of time 
delay on correlated momentum, where the phase angle of 
momentum vector (&[', fejj) is introduced to describe the 
correlation between two electrons' momentum: [0,7r/2] 
and [7r,37r/2] correspond to a same-hemisphere emission 
while [7r/2,7r] and [37r/2,27r] correspond to a opposite- 
hemisphere emission. In all three cases we observe a 
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FIG. 8: (color online). The relationship between the cor- 
related momentum and the delay time at (a)0.12PW/cm 2 
(b)0.4PW/cm 2 (c)1.0PW/cm 2 



FIG. 9: (color online). DI yield versus the delayed time at 
laser intensity of 0.15PW/cm 2 . The total optical cycles is 
(a)7, (b)3, respectively. 



long-tail up to several optical cycles. For the sequen- 
tial ionization of 1 PW/cm 2 , it means after one electron 
is deprived from nuclei by laser fields the other electron 
is slowly (i.e. waiting for up to a few optical cycles) ion- 
ized. In the tunneling regime, the long-tail indicates that 
CEI mechanism is very pronounced for the molecular DI 
(contribute to 80% of total DI yield) . This observation is 
different from purely classical simulation [271 ]. where CI 
effect is believed to be overestimated. Our results, how- 
ever, are consistent with experiments of Ar atom [3l| . 
where the ionization energy and laser field parameters 
are close to our case. The reason is stated as follows. 
For the intensity of 0.12PW/cm 2 , the maximal kinetic 
energy of the returned electron is 3.17U P = 0. 85a.M., still 
smaller than the ionization energy of Nj~ . Even with the 
assistance of the Coulomb focusing [1,13], it is still not 
easy for the returned electrons to induce too many CI 
events. 

Furthermore, such time delay might provide more 
physics beyond simple rescattering scenario. Recently 
a statistical thermalization model has been proposed for 
the nonsequential multiple ionization of atoms in the tun- 
neling regime [33| . This model shows that sharing of ex- 
cess energy between the tunnelled electron and the bound 
electrons takes some time, resulting in a time delay on 
attosecond time scale between recollision and ionization. 
Our simulation upholds this picture of attosecond elec- 
tron thermalization: on upper panel of Fig. [5J two bright 
spots are observed at a similar time delay on subfemtosec- 
ond time scale for CI trajectory. 

The regular pattern in Fig[5Ja),(c) exhibit that 
the ejection of electrons in the same-hemisphere and 
opposite-hemisphere emerge alternately with respect to 
the delayed time. For a time delay of odd half laser cy- 
cles, two electrons are emitted in the same direction, for 



a time delay of even half laser cycles, two electrons are 
emitted back-to-back. Moreover, in the tunneling regime 
the pattern in Fig[5Ja) shows two notably bright spot in 
the first and third quadrants when the delayed time is less 
than 0.5T, a phenomenon directly due to the CI trajec- 
tories. On the other hand, the irregular pattern emerges 
in FigJSJb) for DI in the plateau regime as the signature 
of complicated multiple-collision trajectories. It implies 
that the trajectories of two electrons entangle with each 
other before DI ionization occurs and the electrons' mo- 
tion might be chaotic [3^ ]. 

Integrating over angle in FigJS] gives total DI yield vs 
the delayed time as illustrated in FigOU where a laser 
field with seven and three optical cycles is used in (a) 
and (b), respectively. In each case, we also choose three 
alignment angles, i.e. 0°, 35° and 90° for comparison. 
As known to all, a DI event with delayed time less than 
0.5T represents CI trajectory while others are all CEI tra- 
jectories. In Figl5]Ja), both CI and CEI are suppressed 
for perpendicular molecules compared with parallel ones. 
However, for the case of 35°, CEI play an more impor- 
tant role than other two cases. As a result, the total 
number of CI and CEI events stays almost unchanged 
and forms the quasi-plateau between [0°,35°] in Figj5^a). 
This plateau can eliminated by using shorter laser pulses 
because the excessive CEI can be effectively suppressed 
as demonstrated in FigJ^b). 



V. CONCLUSIONS 

In summary, we have developed a semiclassical quasi- 
static model with including classical rescattering and 
quantum tunneling effects and achieved insight into the 
dynamics of two highly-correlated valence electrons un- 
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der the combined influence of a two-center Coulomb po- 
tential and an intense laser field. Our model can be 
used to compared with experimental data of molecular 
DI quantitatively under the relevant experimental condi- 
tions, i.e., highly nonpcrturbativc fields with femtosecond 
or shorter time resolution. We reproduce the experimen- 
tal data of the plateau structure of the DI, uneven DI 
yield for parallel or perpendicular aligned molecules, and 
the accumulation of the emitted electrons in the first and 
third quadrants of parallel momentum plane. In addi- 
tion, our model calculation reveals some novel features, 
such as a shallow dip in the plateau regime, enhancement 
of anisotropy when illuminated by shorter laser pulses 
and the regular patterns of correlated momentum with 



respect to the time delay. We hope our theory will stim- 
ulate the experimental works in the direction. 
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